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Systematic theoretical results for the effects of a dilute concentration of magnetic impurities on 
the thermodynamic and transport properties in the region around the quantum critical point of a 
ferromagnetic transition are obtained. In the quasi-classical regime, the dynamical spin fluctuations 
enhance the Kondo temperature. This energy scale decreases rapidly in the quantum fluctuation 
regime, where the properties are those of a line of critical points of the multichannel Kondo problem 
with the number of channels increasing as the critical point is approached, except at unattainably 
low temperatures where a single channel wins out. 
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The asymptotic low temperature singularities in the 
thermodynamic and transport properties of many solids 
appear to be due to impurities . While several impurity 
models have quantum-critical points (QCP), where such 
singularities may be expected, they require special sym- 
metries unlikely to be present in real systems An al- 
ternate possibility is that the pure system is near a QCP, 
so that singular low energy fluctuations are present ||-|| . 
We investigate the possibility in this paper that such fluc- 
tuations drive the effects of a dilute concentration of or- 
dinary impurities (not requiring any special symmetries) 
so that the resulting observable properties are much more 
singular than those near pure QCP 



For coupling to non-magnetic impurities, a general 
argument on the renormalized impurity scattering was 
given in relation to insulating behavior in marginal Fermi 
liquids ||. It was recently proposed that critical fluctu- 
ations enhance potential scattering from non-magnetic 
impurities so that the residual resistivity increases as the 
QCP is approached ||. 

Here we investigate the effects of magnetic impurities 
near a ferromagnetic QCP. The problem is especially in- 
teresting for a number of reasons. First, it couples the 
quantum fluctuation or Kondo effect of magnetic impu- 
rities to the singular quantum fluctuations of the pure 
system. Second, since the transition is at q = and the 
order parameter is conserved, conservation laws or Ward 
identities can be used to obtain systematic results for the 
effect of impurities just as they are available for the pure 
ferromagnetic QCP. Third, experimental results on ex- 
traordinary pure samples of a ferromagnet-MnSi M,0, 
which has a QCP as a function of pressure, are in good 
accord with the theory pj. A few results on less pure 
MnSi are available pL which show interesting devia- 
tions from that of the purer samples near the QCP. We 
hope our results will serve as an impetus to further ex- 
perimental results. 



Fig.l presents a schematic phase diagram around a 
pure ferromagnetic QCP at r=0 where r is a "disorder- 
ing" parameter, pressure for MnSi. The crossover be- 
tween the regimes I and II occurs at Toc£ -3 (T) where £ 
is the magnetic correlation length 0. In the pure limit, in 
the quantum-critical regime, the pole in the fluctuation 
spectra %(q, uj) has a dispersion uj^q 3 , so that the dy- 
namical critical exponent Zd—3. Correspondingly, £(T)~ 
j 1 — 2/3 j n -(.jjg re gj on i j n the region II, £~r -1 / 2 at low T. 

The problem of magnetic impurities was looked at long 
ago by Larkin and Mel'nikov |L2|. We have obtained 
some new results. In the regime I, the net result of en- 
hanced dynamical spin fluctuations at low energies and 
of a decreased high frequency cutoff is to enhance the 
Kondo temperature. In the regime II, this enhancement 
decreases rapidly and an ordinary 5=1/2 magnetic im- 
purity necessarily scatters in many angular momentum 
channels because of the increasing magnetic correlation 
length with the number of channels diverging as r— >0. 
The crossover from the multichannel to the single-chan- 
nel behavior occurs at unattainably low temperatures. 

We consider the model of 5=1/2 magnetic impurities 
coupling to the host electrons by the Hamiltonian 
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We assume that the dynamical spin susceptibility of the 
pure system is given by 



x(q»w) = Xo^o^ 2 + q 2 - iw/ft?]" 



(2) 



where xo is the magnetic susceptibility of the noninter- 
acting system, kq is of the order of the Fermi wave num- 
ber pf, and k=k(T) is the inverse of £(T) Q. Hereafter, 
r is defined by r=K 2 (T=0)/4p|. 

For coupling to magnetic impurities, the vertex renor- 
malization AS a \p,vj) / is given by a formally exact 
expression in terms of the irreducible part A^{p, w): 
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where a is the coupling constant between the host elec- 
trons, a=2Uxo f° r the Hubbard interaction. We have 
used the simplified notation, p=(p, s) and zu—(q,u>). 

For fluctuations in which oj~q 3 , a version of Migdal 
theorem holds in the so-called q limit, uj^O, q— >0 with 
uj/vq = so that the Ward- identity requires p| 



A W /A^=[z(m7m) (a + Xo/x)]" 
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i.e., z(m*/m)~a!~ 1 '--'l when Xo/x^l- By general con- 
siderations, m/m*=z(l+d'E/de p ), So <9£/<9e p is not sin- 
gular while z~ 1 = l+dY,/de~\ ln(|e|, T)\. Another conse- 
quence of Eq.(Q) is that the imaginary part of the one- 
interaction irreducible part of xi^) is given by 
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ImxM [fM - /(£k)] S(lu - £ k , + e k )z 2 , (5) 

k 

where k'=k+q, £k is the energy of the quasi-particles. 

The interaction vertex between the host electrons and 
the pseudofermion field has the form r=r(°)+rto -.S. 
To the leading order in J, the vertex is renormalized as 



r (<T) (p,wi;p-|-n7,W2) = 



J AW(p,n7) 
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where lui— l)2=uj. Using Eqs.© and @), is propor- 
tional to x{ m ) m the q limit. In the three-dimensional 
system with a spherical Fermi surface, it is useful to con- 
sider partial- wave components of rto : 



r (<T) (p,£,u/i;p + q,£ + £J,£J2) = 
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ir, (<r) (e,wi;e + w,£^)f}(cos0), 
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where q 2 =(p— p') 2 =p 2 +p' 2 — 2pp'cos 9, Pi are Legendre 
functions of the first kind. The real part of Ij (cj)= 
r^(e,u)i; e-\-u),u)i) has a singular contribution, while 
the imaginary part of rf CT ' (u;) vanishes when lu^O. By 
use of three roots, x\,x% and £3, of a cubic equation, 
x 3 — 2x 2 +x— (cj// n K 3 ) 2 =0, we obtain 
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2N pp' [(x 3 - xi){xi - x 2 ) 

+ 2 cyclic permutations 
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where Ci=(p 2 +p' 2 +K 2 a;i)/2pp' f° r £=1,2,3 and Q; are 
Legendre functions of the second kind. Noting that Qi{() 
have branch points at ^=±1, we can evaluate the singular 
contribution in as 



In 
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where Di=r(2pp exp(— X)L=i «~ 1 )) 3 serves as an upper 
cutoff of \ui\, rn 3 or r|p — p'\ 3 for each Z. For the Kondo 
problem we will investigate bellow, the on-shell scatter- 
ing £p'=£ p +u; is important so that |u;|>_r|p— p'| 3 . Hence, 
we put p=p'=pp. 

If \u)\ is regarded as ~T, Eq.([l0|) leads us to the follow- 
ing facts: in the regime II, the w-dependence of is 
negligible and the number of effective angular momentum 
channels increases as Z max ~2pF£ because the effective size 
of the impurity scales up as £ jL2j; in the regime I, the 
w-dependence of Ij is singular as ocln(.Dj/|w|) be- 
cause the effective exchange interaction acts over the long 
range not only in space but also in time. 

In order to demonstrate carefully how the singular in- 
dependence is combined with perturbative renormaliza- 
tion group (RG) equations for the Kondo problem near 
the QCP, first we consider the renormalization factor of 
the pseudofermion propagator Z(uj). The leading order 
correction, 5Z(lo), with respect to J is given by 



8Z(u) = -j^S{S ■ 
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Noting that the imaginary part of x is given exactly by 



Imx(ro) 
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and using Eqs.(^)-(^), Eq.(pd|) can be evaluated as 
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s(s+i)f: 2 -^ f D ds 
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where g—NmpY /2ir 2 . 

From Eq.(p^), the variation of Z(u>i) on reducing the 
bandwidth cutoff from D to D—6D is obtained immedi- 
ately. For r[ CT '(£ p , uji] e p > , L02) , the 0(J 2 ) term is 
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The 0{ J 3 ) term of TjT' can be evaluated in a similar way 
to 5Z(u). In Eqs.dia) and (0), z and m* should be un- 
derstood as functions of the energy scale D, respectively. 
The high-energy spectrum around the band edge of the 
host electrons may be approximated to be that of the free 
electrons so that z~l and m*^m at an early stage of the 
RG procedure. Even at the low energy, where z vanishes 
on a logarithmic scale, the product z(m* /m) remains of 
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the order of 1 by virtue of Eq. (||) . Therefore the renor- 
malization of the wave function and of the mass of the 
host electron does not affect the resulting scaling equa- 
tion. Neglecting e p /D, e p '/D and lu\/D compared to 
1 while keeping the explicit dependence on u=e p / — e p , 
we obtain the 2-loop scaling equation for the invariant 
coupling Xi (uj) as 



dlnT> 



21' 



(15) 



Z'=0 



Using Eq.(|10|), the initial bandwidth cutoff in Eq.(|lf 
can be put as Z?o=T(2pf) 3 ; the bare coupling con- 
stant, i.e., the value of Xi(lo) at D=D$, is given by 
A| b ''(a;)=(£)/a)r[ <7 ''(a;) and the scaling stops at D=T. 

The solution of Eq.@ has the form, Xi(lj)^X { ^ ] (u)A 
+(Jg/N)Bi, where A and Bi are functions of D which 
are independent of u). Correspondingly, the form of the 
effective time-dependent interaction in the momentum 
space divided by Jg/N, which is given by 



AS e 
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is invariant under the present RG transformation with 
scaling time as t~l/£), where <5o=Ko/4pp and Slg= 
/Jsin((9/2)[sin 2 (6>/2)+<5] with 5=K 2 /4p|. The RG flow of 
A and Bi describes correlation not only between angu- 
lar momentum channels but also between long-time and 
instantaneous components of the effective interaction. 

If we consider the scaling equation at the one-loop level 
where the third-order term of Eq.(|l5|) is neglected, we can 
make a rough estimation of the Kondo temperature Tk, 
which is associated with the breakdown of the perturba- 
tion theory, near the QCP. By use of Eq.([To|) for xj b \uj), 
the one-loop scaling equation can be solved analytically. 
The result depends on whether Tk 3 (T) is smaller than 
T or not in accordance with the region I or II as 
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where J=(J q/6N)(kI/pI). Since A o (0)>A ( (0) for arbi- 
trary i, we evaluate Tk as temperature at which Ao(0) 
diverges. Noting r^c^D^r 3 / 2 in the region II, we obtain 



7 K | exp (-Tr/2y/~jj 
Do I r 3 / 2 exp (-l/v^Jtan 



J|lnr 3 / 2 | 



r < r 
r > r* 



where r*=exp n/Sy/j). As a direct consequence of 
the coupling between the quantum fluctuation of an in- 
dividual magnetic impurity and dynamical spin fluctua- 
tions of host electrons, Tk is seen to be enhanced in the 
region around a ferromagnetic QCP. 



On the basis of the two-loop expansion, the quasi- 
particle's damping rate (r, T)=zlmSi mp (r, T) and 
the electrical resistivity p lmp (r, T) due to magnetic im- 
purities can be expressed in terms of A; (a;) as follows: 

7^ = -S(S + 1) £(i + 1) 2 (A 2 (0) - A 2 +1 (0)), (17) 

1=0 

2 00 

Pimp = ^S(S + 1) J2Q + 1)(^(0) - A i+1 (0)) 2 . (18) 



1=0 



By solving Eq.(jl5|) numerically by use of Eq.® for 



A 



(b) 



j (uj), T-dependences of T imp and p- lrap are obtained in 
Fig. 2, for J=0.02 in the weak coupling regime and for sev- 
eral values of r. We also show the result of the Born ap- 
proximation where Tj^poc^T) and /5 im pCxln£(T), which 
visualizes the crossover line from the regime I to the 
regime II. The fixed point is that of the single-channel 
Kondo problem; the multichannel effects can be seen as 
a transient phenomenon. Note that Tk, at which p lmp is 
of the order of 1, increases rapidly in the region II from 
r~0.1 to ~0.01 while showing a tendency to be saturated 
in the region I, consistently with the one-loop result. 

In the regime II far away from the QCP, the in- 
dependence of the vertex is suppressed, so that 
the low-energy effective Hamiltonian may be mapped 
into the anisotropic n«(/ max +l) 2 -channel Kondo model 
with the band- width cutoff D«D r 3 / 2 , in which the cou- 
pling constant A; is given by A; (0). By general consid- 
erations of the anisotropic multichannel Kondo model, 
while one enters the influence of the multichannel fixed 
point for T below O(Tk), channel anisotropy introduces 
another energy scale T x below which the stable single- 
channel fixed point is reached Hence the crossover 
from the multichannel to the single-channel behavior 
is predicted when T x /Tk<C1. For the n-channel case, 
Tk^DX 71 / 2 exp(— 1/A) and T x asymptotically scales as 
(AA) 1+ ™/ 2 , where A and AA are the average and dif- 
ference of coupling constants, respectively (m). For the 
two-channel case, a recent numerical RG study has shown 
that T X ~T>(AA/A) 2 exp(- 1/A) H). From these results, 
we expect that T x is given by D(AA/A) 1+n/2 exp(-l/A), 
so that T X /T K can be evaluated as (AA) 1+n / 2 /A 1+ ™. 
For the present effective Hamiltonian, from Eq. (fl0|) , 
A;~3J|ln^/r| while A;-A;+i~3 J/(l+l). Since Z max « 
r -1 / 2 , we may make a rough estimation of T x /Tk as 



T 
T K 



- ~ (Sir^lln^-"- 1 a (V3j|lnVF|) 



-l/r 



(19) 



At r=0.1 this is <3(1CT 3 ) and at r=0.05 of 0(1(T 8 ) for 
,7=1. Therefore there is a possibility of T x <§;Tk in the 
intermediate coupling regime. Then for the region of in- 
terest away from the QCP in the region II, at experimen- 
tally attainable temperatures the observable properties 
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are expected to be those for a line of critical points with 
singular properties with exponents continuously chang- 
ing as r changes. 

The true fixed point at r>0 would be the fixed point 
of the single-channel Kondo problem. It has zero ground 
state entropy. For finite r and T^>T X in the intermediate 
coupling regime, we may read off the results for the en- 
tropy and its leading temperature dependence from the 
exact solution for the multichannel problem fl6| . 

On the experimental side, systematic results for the 
resistivity as a function of r are not available. Further 
work is required to test the predictions made above. 
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FIG. 1. A schematic phase diagram near a Ferromag- 
netic Quantum Critical Point. 
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FIG. 2. 



T dependences of r imp and pi mp 



for J=0.02, 



at different r's (r=10" 4 ,10" 3 ,10 -2 , and 1CT 1 going down, 
starting from the top curves at high T) on the basis of the 
two-loop expansion (solid lines) and the Born approxima- 
tion (dotted lines). 
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